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The shock structure of a plasmoid in magnetic reconnection in low-beta plasmas is 
investigated by two-dimensional magnetohydrodynamic simulations. Using a high- 
accuracy code with unprecedented resolution, shocks, discontinuities, and their inter¬ 
sections are resolved and clarified. Contact discontinuities emanate from triple-shock 
intersection points, separating fluids of different origins. Shock-diamonds inside the 
plasmoid appear to decelerate a supersonic flow. New shock-diamonds and a slow 
expansion fan are found inside the Petschek outflow. A sufficient condition for the 
new shock-diamonds and the relevance to astrophysical jets are discussed. 
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I. INTRODUCTION 

Magnetic reconnection provides fast plasma transport in solar, space, laboratory, and 
astrophysical plasmas. By breaking a magnetic topology, it releases the stored magnetic 
energy into the energy of a fast ontflow jet, whose speed is approximated by the Alfven 
speed in the npstream region. The reconnection system often involves a “plasmoid,” a dense 
plasma island conhned in a magnetic loop. For example, when the reconnection jet sweeps 
plasmas and the reconnected magnetic fields, a plasmoid is formed at the jet front. For 
example, the tearing-type instabilities generate plasmoids in a stationary or reconnecting 
current sheet. 

In magnetotail physics, motivated by satellite observations,*^^^ the structure and dy¬ 
namics of a plasmoid have long been studied since the early days of magnetohydrodynamic 
(MHD) simulations. The plasmoid has been considered as an outcome of Petschek-type fast 
reconnection. Carrying out a series of MHD simulations, Ugai^^ extensively studied the 
dynamics of a plasmoid in the downstream of the reconnection outflow. It was found that 
the plasmoid moves at some fraction of the Alfven speed and that it exhibits a complicated 
structure, surrounded by slow-shocks. Abe & Hoshino^ carried out a large-scale MHD simu¬ 
lation and they found a pair of intermediate shocks inside the plasmoid, where the plasmoid 
swallows the thin plasma sheet. Recently, Zenitani &: Miyoshil^ (hereafter referred to as 
Paper I) have introduced a shock-capturing code to the numerical research on the plasmoid 
structure. They found a new class of shocks near the plasmoid. The shocks are attributed 
to the low-beta condition of /? < 1, because the typical speed of the plasmoid motion, the 
Alfven speed (ca), exceeds the local sound speed (c^). These new features suggest that 
plasmoid physics may be qualitatively different in this regime. Since the plasma beta is 
one of the basic parameters in MHD, plasmoid physics in low-beta plasmas deserves further 
investigation. 

Plasmoids draw recent attention in broader context of MHD reconnection. For example, 
Loureiro et al.^ISl predicted that a sufficiently long Sweet-Parker (SP) current sheet is un¬ 
stable to a plasmoid formation. As a result, the current sheet is filled with many plasmoids, 
and then the SP reconnection switches to the “plasmoid-dominated reconnection”®2II5l32^ 
whose reconnection rate is higher than in the SP reconnection. In these studies, the plasma 
beta is usually set to high. However, the plasma beta is extremely low in possible recon- 
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nection sites such as the solar coronaP In low-beta plasmas, the number of the research on 
the plasmoid-dominated reconnection is quite limited, 3^21 and it is totally unclear how the 
above low-beta effects alter the overall evolution. It is of importance to understand plasmoid 
physics in low-beta plasmas as a basic element, in order to explore complex systems such 
as the plasmoid-dominated reconnection. We also note that recent researches suggest that 
the plasmoid-dominated reconnection contains Petschek-like mini-reconnection.l^*^ In such 
a case, the above study is of direct relevance to understand the system. 

The purpose of this paper is to extend our recent work (Paper I) with a better MHD solver 
and to investigate the plasmoid structure in low-beta plasmas in unprecedented detail. We 
discuss several advanced structures in the plasmoid system, such as triple-shock intersections 
and shock-diamonds. 

II. NUMERICAL SETUP 

We use a Godunov-type code to solve resistive MHD equations. The code is based on 
one in Paper I and contains several improvements since then. It employs the HLLD solver.!^ 
Compared with the HLL solver in Paper I, the HLLD solver additionally resolves contact 
and rotational discontinuities inside the Riemann fan. The surface values are interpolated 
by the MUSCL scheme. The minmod limiter is applied to primitive variables. 

The initial setup is similar to the model in Paper 1. The simulations are carried out in 
the x-z two-dimensional plane. We set the reconnection point to (x, z) = (0,0). We solve 
one quadrant of [0, 200] x [0,150] with 12000 x 9000 grid points. The resolution is twice 
better than in Paper 1. We use a symmetric boundary at x = 200 and a reflecting wall at 
z = 150. We consider a Harris sheet: B{z) = tanh( 2 :)*, n = 0, p{z) = 1 -f cosh“^( 2 ;)//dup, 
p{z) = 0.5(/3up + cosh“^( 2 :)), where /Sup is the plasma beta in the upstream region. These 
symbols have their standard meanings, unless stated otherwise. The Alfven velocity in the 
upstream region is set to CA,up = \B\/^/f> = 1. Thus the unit time corresponds to the Alfven 
transit time across the Harris sheet. We set the polytropic index to 7 = 5/3. 

We employ a localized resistivity, ^(x, z) = po + (71 — po) cosh“^(\/x^ -|- where 71 = 
1/60 and 70 = 1/1000. A weak perturbation 5Ay = —0.06exp[—(x^ -h z‘^)/4] is imposed 
to magnetic helds to help the reconnection onset. In this way, we invoke Petschek-type 
reconnectiorP^ to obtain a single, giant plasmoid at the jet front. We intend to understand 
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basic structures in this idealized configuration. We study two representative cases: /dup = 0.1 
in Run 1 and /Sup = 0.2 in Run 2. 

III. RESULTS 

We examine Run 2, which corresponds to the main run in Paper I. Figure is a snapshot 
at t = 250. The upper half shows the outflow velocity and the lower half shows the 
divergence of the velocity field Vv. As discussed in Paper I, there are characteristic shocks, 
such as a Petschek slow shock (PK), a reverse fast shock (FS), a postplasmoid vertical slow 
shock (SS-1), and a forward vertical slow shock (SS-2). In the lower half, these shocks 
are evident in compressional layers of V-t; < 0. The postplasmoid vertical slow shock is 
sometimes found on the backside of the plasmoid. Such configuration is similar to a normal 
shock on the airfoil, when an aircraft is flying at a transonic speed.^ These shocks are features 
of low-beta reconnection, in which an Alfvenic outflow drives the plasmoid at a transonic or 
supersonic speed. 

The forward slow-shock connects with the plasmoid boundary. Since the plasmoid bound¬ 
ary is a slow shock,^ it is interesting to see how three shocks intersect each other. Shown 
in Figure [^a) is the plasma density near the forward slow-shock. The dense (orange) re¬ 
gion corresponds to the plasmoid. The forward slow-shock meets the plasmoid boundary at 
{x,z) ^ (102,8.5). 

In hydrodynamics, we expect another shock or discontinuity from such a triple-shock 
intersection.!^ Consider the fluid advection from the upstream A region to the downstream 
plasmoid. There is one shock between in A —)■ C", while there are two in A —)■ U. Since 
the Rankine-Hugoniot relation is nonlinear, the fluid states should be different in C and 
C. Then a shock or a discontinuity should separate the two regions. We expect that this 
argument holds true in MHD, because MHD is a complex extension of hydrodynamics. In 
this case, there is a contact discontinuity (CD), as indicated by the white dashed line in 
Figure l^a). The specific entropy {p /not shown) is useful to to visualize it. 

Shown in Figure |^b) is the plasma density in the outflow region behind the plasmoid. 
The fast shock is located at x ~ 61. The plasma density decreases in the upstream side of 
the shock, because an adiabatic-type acceleration takes place there.!23 As illustrated, there 
are two triple-shock intersections. We similarly expect additional shocks or discontinuities. 
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A contact discontinuity stems from the left intersection point, where the vertical slow-shock 
connects with the Petschek shock. This CD was suggested, but unclear in Paper I. In this 
work, the HLLD solver successfully resolves it between the C and C’ regions. The CD 
survives even in the downstream of the fast shock, as illustrated by the lower dotted line 
in the diagram. In addition, we hnd another CD from the right intersection point toward 
the upper-left direction. As a result, between the two discontinuities, a narrow low-density 
channel is formed inside the plasmoid. The divergence V v does not mark them in the lower 
half of Figure This suggests that these structures are discontinuities rather than shocks. 
These discontinuities inside the plasmoid have never been discussed in previous works. 

We also hnd a very weak compressional layer at x = 105, as indicated by the small 
arrow in Figure [Tj It sits in the front side of the normal shock since the earlier stage of the 
simulation. Unfortunately the layer is too faint to obtain shock quantities. For example, 
our Rankine-Hugoniot analysis tells us that slow-mode Mach numbers are around 2 in both 
upstream (right) and downstream (left) regions, but it is difficult to evaluate the shock- 
propagating speed. At this point, we are not sure whether it is a compressional layer that 
could evolve into a slow shock, or the HLLD solver does not adequately resolve a slow shock. 

Next, we discuss another key player, shock-diamonds. The plasmoid exhibits a crab-claw 
structure, because a dense plasma sheet lies at the center. When the bifurcated flows are 
faster than the sound speed inside the plasma sheet, the twin edges invoke oblique shocks, 
resulting in a diamond-pattern inside the plasma sheet.®® Such shock-diamonds are evident 
in the right side in Figure (see also Fig. 5b in Paper I). 

In the rest frame of the plasmoid, a narrow plasma sheet comes in at a supersonic speed. 
This is equivalent to a supersonic nozzle problem in aerospace engineering.*^*^ Let us consider 
that a supersonic flow emanates from a nozzle, as illustrated in Figure It is classihed to 
over-expanded and under-expanded cases. If the jet pressure is lower than the ambient 
pressure, the flow is said to be over-expanded. In this case, the oblique shock propagates 
inward from the nozzle. The jet boundary reflects the oblique shock as an expansion wave, 
and then reflects the expansion wave as a compressional (shock) wave. This cycle results 
in a diamond-shaped pattern in the downstream. When the jet is under-expanded, i.e., at 
higher pressure than the ambient pressure, an expansion wave appears hrst, and then the 
same cycle continues. The bottom panel in Figure shows the flow pattern in this cycle. 
The flow converges in the downstream of the oblique shock, and diverges after crossing 
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the high-pressure cell. Due to the pressure difference, the jet boundary also exhibits cyclic 
expansion/compression. The sign of the vertical velocity changes accordingly, as indicated 
by blue and red arrows. 

Next we examine the role of the shock-diamonds. Figure presents ID cuts of selected 
quantities at the midplane, z = 0. They exhibits periodic oscillation in the shock-diamond 
region of 95 < x < 120. The magenta curve in Figure [^indicates the sound speed toward the 
+x direction. Note that the sound speed is virtually identical to the fast-mode speed in the 
right side of the tangential discontinuity (TD) at x ~ 80, which separates the reconnected 
magnetic held and the initial plasma sheet. It is Cg ~ 0.42 in the rightmost plasma sheet. 
This is also indicated by the dashed line in magenta for comparison. The propagation speed 
of the TD is utd ~ 0.54, similar to the plasma speed around the TD. Since it exceeds 0.42, 
macroscopically, the plasmoid interacts with the plasma sheet at a supersonic speed. In the 
plateau of 85 < x < 95 the velocity is Ua, ~ 0.1 and so the right-running sound speed is 
{vx + Cg) ~ 0.58 > utd- Between the plateau and the leftmost plasma-sheet, we recognize 
gradual slopes in how properties, as indicated by the dashed lines. These results suggest 
that the shock-train decelerates the plasma-sheet how from a supersonic speed to a subsonic 
speed in the comoving frame of the plasmoid. 

In huid dynamics, it is known that such a shock-train decelerates a supersonic how.h^*^ 
Similar examples include a supersonic Fa n n o how and an oblique shock dihuser. In Figure 

the central shock channel becomes narrower and the Mach angle of each shocklet becomes 
steeper, as the plasma sheet penetrates deeper into the plasmoid. The typical length of 
the shock-train region is 10-20 times long as the how width.^ These are common features 
of the shock-train that decelerates a supersonic how. From the right upstream to the left 
downstream in Figure]^ the pressure increases in the shock-diamond region. This is mainly 
due to adiabatic heating. The strongest shock in the shock-train is the rightmost oblique 
shock, but its Mach number is only A4 ~ 1.2. The entropy gain is limited across such a weak 
hydrodynamic shock. It appears that the shock-train works as a quasi-adiabatic decelerator. 

We expect shock-diamonds, wherever a narrow jet crosses a boundary at a supersonic 
speed. Motivated by this, we hnd another shock-diamonds. Figure |^a) shows the vertical 
velocity Vz near the intersection of the Petschek shock and the normal shock at t = 275 in 
Run 1 (see also Fig. 6(a) in Paper I). Note that the range is carefully set to —0.02 < Vz < 0.02 
in the unit of CA,up- The normal shock is indicated by the white dash line. As evident in 
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the triangle pattern in one can recognize new shock-diamonds inside the reconnection 
jet. They are of under-expanded type, because the pressure in the B region is higher than 
in the A region, i.e., > Pa- As color-coded in Figure |^a) and as illustrated in Figure 

[^a’), the upward motion in the expansion phase comes hrst in Figure]^ (see also Fig. [^). 
The vertical speed is small, \vz\ ~ 0.01, and the fluctuation in pressure is Ap/p < 0.1. The 
shock-diamonds are marginally visible in Run 2, but they are much more apparent in Run 
1. They further evolve longer in larger simulations. 

Let us examine the condition for the new shock-diamonds. The shock-diamonds will 
appear when the reconnection jet supersonically outruns the normal shock, Vjet ~ CA,up > 
Csjet + "yshock, where ca is the Alfven speed and Cs is the sound speed and the subscripts 
“up,” “jet,” and “shock” denote the upstream region, the outflow jet, and the normal shock, 
respectively. We note that plasmas in the B region are stationary and its properties are 
similar to the initial upstream plasmas. Since the normal shock is a slow shock in a low-beta 
plasma, we expect that the shock-propagating speed is subsonic, Vshock < Cg^up- Then the 
sufficient condition for shock-diamonds is 


C-AjUp ^ \ i T 


Pjet V dup 

We further utilize the pressure balance across the Petschek shock, 

1 + /3up 


( 1 ) 


Pjet ~Pup- 


A 


+ Pup(0.1CA, 


2 

upy 1 


( 2 ) 


up 


where the last term is negligible, ~ (pup//9up)Cl(10 ^). The compression ratio is known to 

btpa 

Pjet _ F(1 -F /3up) 


Pup 1 + F /3up 

where F = 7/(7 — 1). Using these relations. Equation ([^ yields 


(3) 


2F 


a/FAup + a/ 1 + FAup <1/ -• 
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This gives a simple condition. 


Aup < 0.135. 


(4) 


(5) 


Run 1, Aup = 0.1, is in this regime. In the Aup ^ 1 limit, shock-diamonds will always appear 
inside the Petscheck outflow. 
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Surprisingly, the post-plasmoid structure is more than a superposition of the contact dis¬ 
continuity (CD) and the shock-diamonds. Figures [^b) and (c) show the plasma outflow 
and the specihc entropy s = p/p'^. Our interpretation is illustrated in Figure [^b’). The en¬ 
tropy, which is an excellent measure of the fluid state, visualize the CD from the triple-shock 
intersection in the horizontal direction. This is indicated by the white dotted line. Interest¬ 
ingly, the held lines are bent below the CD. The bent point propagates inward, as indicated 
in the dot-dash line. We think that this is a slow expansion wave (SE). The C" region 
between the CD and the wavefront corresponds to the slow expansion fan. The following 
results further support our interpretation. The narrow angle of the C" region tells us that 
the wavefront propagates at the local Alfven speed or slow-mode speed. Possible candidates 
are slow shock, slow expansion wave, and rotational discontinuity. From the upstream (C") 
to the downstream (C"'), the magnetic tilt angle changes from ?«90° to ?«50°. Since the 
wavefront angle is only about —5°, the rotational discontinuity is ruled out. Since both the 
tangential magnetic held and the magnetic held strength increase across the wavefront, the 
slow shock is ruled out. Therefore we conclude that this is a switch-on type slow expansion 
wave. Since plasma beta is high in the outhow exhaust, changes in plasma properties are 
rather subtle to distinguish. The entropy remains unchanged (Fig. |^c)) in C and in C'\ as 
it should be. 


The formation mechanism of the expansion fan can be understood as follows. In the 
upstream side of the Petschek shock, Alfven speeds are ca = 1.05-1.08 in the rarehed 
A region and ca = 0.9 in the shocked B region. Accordingly, from the Petschek balance 
condition, we expect that the outflow speed in C’ is slightly faster than in C. These flows are 
separated by the CD. In hydrodynamics, the tangential velocities can be different across the 
CD. This is why the CD is often called the “slip line” in the context of the shock theory.l^ 
However, in this MHD case, the two flows are threaded by the same magnetic held lines 
inside the outhow region. So, as evident in Vx in Fig. |^b), the fast how in C’ drags plasmas 
in C” via the magnetic held line. This explains why the held lines are bent at the wavefront 
and why the C” region is slightly rarehed. 
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IV. DISCUSSION 

Thanks to the HLLD solver and the very high resolntion, we have identified contact 
discontinnities near the triple-shock intersection points. We have also examined two shock- 
diamonds. One is inside the plasmoid and works as an adiabatic decelerator for the incoming 
plasma-sheet flow. The other is newly fonnd inside the Petschek ontfiow. The post-plasmoid 
strnctnre fnrther involves the slow expansion fan. Many of them arises from intersections 
between the classical Petschek slow shocks and the normal shocks, the latter of which appear 
in the (3 <1 regime (Paper I). In MHD, the discontinnities and slow-mode strnctnres develop 
and persists, nnless they are blnrred by the nonideal effects. At this stage of investigation, 
their impact on the overall evolntion is not clear, however, the system certainly becomes 
more complex than expected. They were not adeqnately resolved by the HLL solver, bnt 
the HLLD solver snccessfnlly resolves them. 

We believe that these featnres are nbiqnitons in low-beta reconnection, as long as the 
MHD approximation holds trne. While this work has focused on Petschek reconnection, the 
plasmoid interacts with the plasma sheet either in Sweet-Parker and plasmoid-dominated 
reconnections as well. Once a plasmoid moves at a sufficiently high speed, it is followed 
by normal shocks in whatever form of reconnection (e.g.. Fig. 3 in Ref. I33|) . Indeed, in our 
preliminary results on plasmoid-dominated reconnection, we confirm that plasmoids invokes 
a lot of normal shocks. This makes the upstream region highly dominated with shocks. 
Furthermore, recent researches report that the plasmoid-dominated reconnection involves 
Petschek-type structure with a pair of slow shocksFinally, during the turbulent regime, 
a plasmoid not only hits the plasma sheet ahead of it, but also collides into another plasmoid. 
All these considerations suggest that shocks and discontinuities further intersect with each 
other in the plasmoid-dominated regime. 

We are aware of some limitations in our model. From the application viewpoint, we have 
to consider several non-MHD effects. For example, the CD will be blurred by collisionless 
mixing along the field lines in the magnetotail. In solar corona, the heat conduction along 
the field lines, radiative cooling, and collision with the neutral atoms will modify the sce¬ 
nario. Our results will be a baseline to discuss these non-MHD effects. From the numerical 
viewpoint, the HLLD solver neglects slow waves inside the Riemann fan. Therefore, it basi¬ 
cally tends to blur slow-mode structures. In particular, in a low-beta plasma, the fast-mode 


9 


PHYSICS OF PLASMAS 


and Alvfen speeds are much higher than the slow-mode speed in the field-aligned direction. 
This is one of the reasons why the compressional layer is ambiguous in Section 3. The slow 
expansion fan is obviously beyond the scope of the HLLD solver. To resolve these slow-shock 
structures, we may need to employ an exact Riemann solver.®^ 

It is interesting to mention that our single plasmoid system is somewhat similar to the 
astrophysical jet system. Figure schematically illustrates a typical structure of the jet.^ 
Aside from the minor difference between the upper and lower panels, we outline the essence 
in the following way. A narrowly-collimated jet travels at a supersonic speed, hits the inter¬ 
stellar matter, form a reverse shock in front of the CD, and then turns round as a backflow 
inside the cocoon. The supersonic jet is thought to contain oblique shocks, which may be 
observed as periodic knots.^^ The oblique shocks are attributed to either shock-diamonds 
or similar periodic structure from the jet-launching sit^SElIZElIlS] shocks driven by the 
jet-cocoon interaction. In our system, in the comoving frame of the plasmoid, the nar¬ 
row plasma-sheet travels at a supersonic speed, hits a dense plasma region in front of the 
tangential discontinuity, and then turns round as a backflow inside the plasmoid. Interest¬ 
ingly, we do not recognize a reverse shock inside the plasmoid. Plasmas are compressed 
around 80 < x < 84 in Figure but there is a gradual slope in front to the TD. This 
is because the flow is already decelerated by the shock-train. The reverse shock does not 
stand, because the TD propagation speed is slower than the the right-running sound speed, 
{vx + Cs) ~ 0.58 > utd- We attribute this to over-expanded shock-diamonds. As the ambi¬ 
ent medium confines the jet, a supersonic flow tends to be adiabatically decelerated in this 
case. This will also enhance the interaction between the jet and the surrounding medium. 
In contrast, the astrophysical jet is usually under-expanded.^ Since released to a tenuous 
space, a supersonic jet can be adiabatically accelerated. It is reasonable that the deceleration 
effect is less prominent. In summary, the plasmoid structure is similar to the astrophysical 
jet structure, but it lacks a reverse shock due to the deceleration effect in overexpanded 
shock-diamonds. 

In this work, we have explored advanced shock-structures around the plasmoid in mag¬ 
netic reconnection in low-beta plasmas, such as triple-shock intersections and two shock- 
diamonds. The new structures will be basic elements to discuss various reconnection- 
plasmoid systems in low-beta plasmas. From the viewpoint of basic fluid physics, these 
structures are the outcome of compressible fluid effects or high-speed fluid dynamics.An 
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adiabatic effedP^ and a new class of shocks (Paper I) also belong to compressible effects. 
These compressible effects appear, when the typical flow speed is comparable with or is 
faster than the sonnd speed. Low-beta plasmas nsually satisfy this condition, becanse the 
Alfvenic motion easily becomes transonic or snpersonic. The compressible flnid dynam¬ 
ics will be a key to nnderstand magnetic reconnection and associated MHD phenomena in 
low-beta plasmas, inclnding the plasmoid-dominated reconnection. 
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FIG. 1. (Color online) Snapshots at t = 250 in Run 2, overlaid by magnetic field lines. (Top) 
Plasma outflow velocity Vx- (Bottom) The divergence of the velocity field V-v, magnified by 5. 
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FIG. 2. (Color online) Plasma density p and schematic drawings in (a) the forward slow-shock 
region, and (b) the postplasmoid region. 
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(a) Over-expanded flow 



(b) Under-expanded flow 



(c) Flow pattern 
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FIG. 3. (Color online) Supersonic nozzle problem, (a) Over-expanded flow, (b) under-expanded 
flow, and (c) the flow pattern in the shock cell. 
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FIG. 4. (Color online) ID cuts of MHD quantities at z = 0. The plasma density (p), the pressure 
(p), the right-running sound speed {vx+Cs; only at x > 80), and the plasma speed (vx) are preseted. 
The ‘TD’ indicates the position of the tangential discontinuity. 
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(b) Outflow velocity, Vx xio° 
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FIG. 5. (Color online) Properties near the shock intersection point in Run 1, (a) vertical velocity 
Vz, (b) outflow velocity Vx, and (c) specific entropy s = p/■ The vertical dash lines indicate the 
normal shocks. The bottom two panels are overlaid by the magnetic field lines. In (a) and (b), the 
ranges are carefully adjusted to emphasize the structure. The right panels schematically illustrate 
(a’) shock-diamonds and (b’) discontinuities. 
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FIG. 6. (Color online) Structure of the astrophysical jet. (Reproduced from Mizuta et al.^; 
Courtesy of A. Mizuta) 
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